Stars and (Furry) Black Holes in Lorentz Breaking Massive Gravity 



D. Comelli, 1 F. Nesti, 2 and L. Pilo 3 

1 INFN, Sezione di Ferrara, 1-35131 Ferrara, Italy 
2 Dipartimento di Fisica, Universita di Ferrara, F44100 Ferrara, Italy 
3 Dipartimento di Fisica, Universita di L'Aquila, 1-67010 L'Aquila, Italy, and 
INFN, Laboratori Nazionali del Gran Sasso, 1-67010 Assergi, Italy 



We study the exact spherically symmetric solutions in a class of Lorentz-breaking massive 
gravity theories, using the effective-theory approach where the graviton mass is generated by 
the interaction with a suitable set of Stiickelberg fields. We find explicitly the exact black-hole 
solutions which generalizes the familiar Schwarzschild one, which shows a nonanalytic hair 
in the form of a powerlike term r 7 . For realistic self-gravitating bodies, we find interesting 
features, linked to the effective violation of the Gauss law: i) the total gravitational mass 
appearing in the standard 1 jr term gets a multiplicative renormalization proportional to the 
area of the body itself; ii) the magnitude of the powerlike hairy correction is also linked to 
size of the body. The novel features can be ascribed to the presence of the Goldstones fluid 
turned on by matter inside the body; its equation of state approaching that of dark energy 
near the center. The Goldstones fluid also changes the matter equilibrium pressure, leading 
to an upper limit for the graviton mass, m < io~ 28H ~ 29 eV, derived from the largest stable 
gravitational bound states in the Universe. 
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I. INTRODUCTION 

The question whether general relativity (GR) is an isolated theory is interesting from both the 
theoretical and phenomenological side. It is known that one can add to the Einstein-Hilbert 
Lagrangian a tower of higher order operators, starting with terms quadratic in the curvature, 
which produces corrections ~ (E/Mp) n , where Mp is the Planck mass and E is the typical energy 
of the process we are considering. Thus the predictions of GR at scales much larger then the 
Planck length 1 /Mp are insensitive to them. Modifying GR at large distances is a totally different 
business. Besides its theoretical interest, the search for large-distance modified theories of gravity 
has been motivated by the evidence for cosmic acceleration and the consequent revival of the long- 
standing cosmological constant problem. The idea is to look for massivelike deformation of GR 
featuring a large-distance (infrared) modification of the Newtonian gravitational potential and a 
massive graviton. 

It is instructive to start from perturbative gravity by considering a Lorentz-invariant theory of 
a massive spin- two field p] . The resulting theory is plagued by a number of diseases that make it 
probably unphysical. First, the modification of the Newtonian potentials is not continuous in the 
limit of very small graviton mass m giving a large correction (25%) to the light deflection from 
the sun that is experimentally excluded [2]. The discontinuity manifests itself in the weak-field 
regime and a possible way to circumvent the problem at the full nonlinear order was proposed 
by Vainshtein in [BJ: if in the Fierz-Pauli theory (FP) the linearized approximation breaks down 
near the Sun the above mentioned discrepancy cannot be trusted anymore. He proposed using an 
improved perturbative expansion, which has a continuous limit for m — > 0. The relative solution is 
valid only up to a finite distance, and the question is then whether this solution can be extended 
up to infinity and matched with the Yukawa-like solution valid at large distances [8]. Recent 
evidences that this can indeed be achieved by means of a different weak coupling expansion have 
been addressed in [3]. 1 In addition, the FP theory is also problematic as an effective theory at 
the quantum level. Regarding FP as a gauge theory where the gauge symmetry is broken by a 
explicit mass term m, one would expect a cutoff A2 ~ mg^ 1 = (mMp) 1 / 2 ; however, the real cutoff 
is A5 = (m 4 Mp) 1//5 [5j much lower than A2. A would-be Goldstone mode is responsible for the 
extreme UV sensitivity of the FP theory, which becomes totally unreliable in the absence of proper 
UV completion. These issues cast a shadow on the possibility of realizing a Lorentz-invariant 
theory of massive gravity [9]. 

It was recently noted that by allowing Lorentz-breaking mass terms for the graviton the result- 
ing theory can be physically viable [10], being free from pathologies such as ghosts or low strong 
coupling scales, and still lead to modified gravity. Since mass terms break the diffeomorphism 
invariance anyway, this possibility was analyzed mainly in a model-independent way, by reintro- 
ducing the Goldstone fields of the broken gauge invariance; and by studying their dynamics [5} HI]. 
For a recent review, see [llj. Lorentz-breaking massive gravity was also considered in the frame- 



1 It is important to anticipate here that also the exact solutions obtained below are free of the discontinuity pr oblem 
and can be obtained [15] by a non-canonical weak coupling expansion where hoi ~ hoo ~ hij ~ e, see Section IIIB 



3 



work of bigravity [12]. In particular, we will be interested in a special phase of Lorentz- breaking 
massive gravity that has no propagating scalar perturbations around flat space. In this phase, the 
absence of dangerous instabilities survive in the curved spacetime, as shown in |13j . 

Besides its consistency, in order to be a viable theory massive gravity has to pass a number a 
tests. In GR the Schwarzschild solution is a benchmark and it is thus crucial to study the impact 
of a massive deformation on it. The main goal of this paper is to study both analytically and 
numerically the gravitational field produced in by a spherically symmetric body. The outline is the 
following. After the quick definition of the theory in the Stueckelberg approach in Sec. |II| in Sec. Ill 



we find the generalized Schwarzschild solution for Lorentz-breaking massive gravity, and explicitly 
obtain the values of the two integration constants entering the exterior solution as a function of 
the mass and the size of the body. This result is obtained matching the exterior with the interior 
solutions for an object of constant density and using the nonstandard perturbative expansion. In 



Sec. IV we discuss the validity of perturbation theory and present a numerical analysis supporting 
our results also when gravity is strong inside the body. In section [V] we discuss some features of 
the black-hole solutions. 



II. IR MODIFIED THEORIES 

Infrared modified gravity theories with manifest diffeomorphism invariance can be realized by 
introducing a set of four Stiickelberg fields & A (A = 1, . . . , 4) transforming under a diff 5x^ = ^(x) 
as simple scalars. Then one can build new geometric objects, manifestly diff invariant, function 
of the metric field and the scalar fields, g^ u d^<& A d v <& B . A generic function of these quantities, 
when expanded around a background, can give rise to graviton mass terms. One can require 
that expanding around a Minkowski background the resulting mass terms are Lorentz-invariant. 
However, as discussed in the introduction, Lorentz invariant massive gravity is rather problematic; 
therefore, we relax this requirement and keep just rotational invariance in the theory. This approach 
is suited to describe a generic theory in an effective fashion, and a general action can be written as 

S = [ d^x^-g (Mp R + £ matt ) + S goki , 

(!) 

S g oid = J d A x^m 2 M 2 P F{X, V\SV) , 
where J 7 is a rotationally invariant potential, function of 

X = -g^d^°d u ^° , V* = -g^dft&dvQ , S ij = -g^d^d^i . (2) 

The constant m sets the graviton mass scale. Note that the breaking of Lorentz symmetry is 
directly built in the action; in Appendix [A] we discuss a different approach where the breaking is 
dynamical. Furthermore, additional symmetries of the Goldstone action can be used to single out 
a particular phase of massive gravity |14j . In particular taking 



T = F{X,W lj ), 



(3) 
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where W % i = S 13 ' — X~ l V 1 V J , the Goldstone action is invariant under <3? 1 — > <& l +£ l ((/> ) and it turns 
out that, in a flat background, only the massive spin 2 tensor (2 elicities only) propagates |14j. 2 
The flat background admitted by the action ([I]) can be parametrized as 



g^u = V^u, $ = (at, bx l ). (4) 

The background breaks boosts and preserves rotations, as can be seen by considering the additional 
effective background metric gi^y = d^ A d u ^ B r]AB = diag(— a 2 , b 2 , b 2 , b 2 ) (see Appendix |a| . The 
condition for the existence of the flat solution Q is the vanishing of the background Goldstones 



energy momentum tensor [see (Bl) in Appendix |B| . Due to rotational symmetry of the potential 
and of the effective g~2n,v > this condition amounts to two independent equations that determine a, 
b in terms of the parameters entering the potential. The background breaks Lorentz when a ^ b, 
but as noted above, the Lorentz breaking is built into the potential, not only in the background 
configuration. This is to be contrasted to theories with spontaneous breaking of Lorentz symmetry, 
as for instance in bigravity. 3 Here, it is notable that for a generic potential function T a flat solution 
is always present, regardless of a cosmological constant term in the action. 

Because of the nonzero background of the scalar fields, their fluctuations with respect to the 
background, tt a = <fr A — Q A ; trasform under diffs as a Goldstone field, 5tt a = ^(x)d^ A /dx^. 
A choice of coordinates setting the Goldstones fields tt a = fixes the gauge completely and is 
sometimes called the unitary gauge: all the dynamics is transferred to the metric. 



III. SPHERICALLY SYMMETRIC SOLUTIONS 

For the spherically symmetric case, one can always find a set of coordinates where the metric and 
the Goldstone fields have the following form 

ds 2 = -dt 2 J(r) + K(r) dr 2 + r 2 dti 2 , 

i (5) 

<S>° = at + h(r), = — 6( r ). 

r 

We will first derive the exact vacuum spherically symmetric solution describing the exterior part 
of the self-gravitating body, then we look for the interior part. Finally, the exterior and interior 
parts of the solution are matched to find the two integration constants, M and 5m which account 
for the gravitational mass of the body and the size of the powerlike term in the solution, one of 
the main novel features compared to GR. 



2 It is quite interesting that a non perturbative realization of this well-behaved "phase" is realized also by a general 
bimetric theory [12] . where Lorentz symmetry is broken by the two tensor condensates, in the limit where the 
second metric is decoupled. 

3 When a = b there are still two equations to be solved, while in bigravity (see appendix [aJ a Lorentz invariant 
background gives only one equation. In this respect the bigravity approach is preferable. 
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A. Exterior, vacuum, exact solution 



The Einstein equations in vacuum read 

K = T e M i (6) 

2M 2 gu ' V ; 

where E„ is the Einstein tensor and T g ^ is the Goldstone energy momentum tensor (see Appendix [b| 
for details). Some general features can be established independently from the choice of T . From 
the ansatz ([5]) it is clear that the Einstein tensor is diagonal and therefore also the Goldstone 
energy momentum must be diagonal: T s t = 0. This gives 

' V [X 4 T X (JKf + X 2 T\4> 2 {JKf - 2XT 2 </»' 4 JK + 3T 3 0' 6 ] = (7) 



K 7 X A 



where T\ ,T 2 ,T 3 and T x are defined in Appendix [Bj and X [see (C6)] depends on h', J ,K. As 
a result, two branches arise: in the first the equation is solved by b! = 0; in the second, with 
h' 7^ 0, the term in square brackets of ^ can be solved for X, or e.g. cp' . As noted already in the 
bigravity context, the first branch h! = 0, leads to equations that are very hard to solve analytically; 
therefore, we will concentrate on the second branch. 

Then, by using the solution to T g * = 0, one remarkably finds that T g * — T g T r = independently 
from the potential. Therefore E\ — E r T = 0, which implies that, as in GR, 

where ko is an integration constant. It not difficult to see that the remaining Einstein equations 
consist in 2 equations for 2 unknown functions J(r), h(r), that can in principle be solved. 

To proceed further an explicit form for T must be provided. A quite general choice of J-, 
inspired to a class of bigravity theories studied in [15] , is 



+ a + a\wi + a 2 (wf - w 2 ) + a 3 (wf - 3w 2 wi + 2w 3 ) . 



X~ x + 



(9) 



where w n = Tv(W n ). The truly remarkable feature of the class of potentials ^ is that the equation 
T g * = admits (still for nonvanishing h') a simple solution with 

(f> = br, ai-4a 2 + 6a 3 - ^o + 2^i -2^ = 0, (10) 



where we defined a n = b~ 2n a n , f3 n = a~ 2 b 2n (3 n . The condition (10) should be regarded as 
an equation determining the Goldstone background a (or b) in terms of the parameters of the 
potential. Using (10) in the remaining equations, one obtains h' 2 (see Appendix [cj and finally 



the exact "black-hole" solution 

2GM 



J(r) = l- — +AV + 2GSr 7 . (11) 
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Here M, S are two integration constants, while G = l/WnMp is the Newton constant and A 2 = 
\m 2 (l2a3 — 60:2 + «o ~~ 3/3i + 12/^2 ~~ 18/33 ) is an effective "cosmological constant". The exponent 
7 is given by 7 = — 2 (2a2 — 60:3 + Pi — 2^2) /{oti — 4q?2 + 60:3). This kind of solution was first 
found in [15] in the context of bigravity theories. 

While the 1/r qnd r 2 terms are also present in the Schwarzschild-de Sitter solution, the last 
one represents a power-law correction to GR whose size is the new integration constant S. It does 
not contain the graviton mass scale, therefore this term is present in the exterior solution also in 
the limit of vanishing graviton mass. In other words the influence of the Goldstone modes survives 
even in the limit m — > 0. 

Clearly, one can choose b such that the effective cosmological constant vanishes, A = 0. Then, 
for 7 < 2, the metric describes an asymptotically flat space, which is just the flat solution Q. The 



condition A = together with (10) are the two conditions for a flat solution determining a, b. The 



masses of fluctuations around flat space (see [10]) are 



O 3 / — — — \ o o q1/ — — \0 

m = - [ax -4a 2 +6a 3 -/3i +4/3 2 — 6/3 3 J m , m x = , m 2 = - (<5i -2a 2 -2/3 2 ) m , 

m| = - (ai -6a 3 -ft +8^2 -18/9 3 ) ^ 2 , m 4 = ^ («i ~ 4 «2 +6a 3 -3^i +12^2 -18^3) m 2 (12) 

and since the mass of the (spin- two) graviton is m|, the corresponding combination of coupling 
constants should be positive. 

For 7 > — 1 the new term S*r 7 is dominant over the Newtonian term M/r, and accordingly the 
total gravitational energy of the solution, evaluated via the Komar mass integral, is infinite (see 
Appendix [E]) . For 7 < — 1 instead the new term is subleading and the Komar energy is just the 
mass M. In the following we will limit ourselves to the case 7 < — 1. 

If the metric represents the exterior geometry of a spherical "star" ; the integration constants M 
and S can be computed in terms of the star parameters by matching the interior and the exterior 
solutions. This will be done in the next section. 

Having the exact solution, it is interesting to discuss its behavior at large distances, in the case of 
asymptotically flat solutions, where it can be compared to the standard weak-field limit. Curiously, 
while J (and K) are of the form 1 + e, with e ~ 1/r, the expression for h', shows that asymptotically 
h ~ y/e. So, the solution does not fit in the standard democratic weak-field expansion, or in other 
words the standard weak-field expansion fails to capture the asymptotic behavior. The correct 
weak-field expansion is actually nondemocratic, where some fields are smaller than others. By 

a gauge transformation one can eliminate h' and turn on the hoi asymptotic components of the 

1/2 1/2 

metric, that again are nondemocratic: hoi ~ h Q0 ~ ; this is analogous to the suggestion of 
Vainshtein for the Fierz-Pauli (Lorentz-invariant) theory, where in the vicinity of a source this 
expansion leads to a non analytic solution which is continuous in the graviton mass parameter 
m — > 0. In that case however, the non analytic solution is valid only up to a finite radius. The 
standard weak-field expansion, leading to the Yukawa falloff is in turn valid only at larger distances, 



4 This was first pointed out in [15] in the context of bigravity. 
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but is discontinuous in m. It was the aim of [3] to setup a nondemocratic expansion to match the 
large and small distance solutions. Here, it is remarkable that the above nondemocratic weak-field 
expansion is valid at all distances, as shown by the exact solution. In other words, the theory is 
always inside the Vainshtein radius. 



B. Inside a -Star- and matching 

Consider now a region where matter is present in the form of a perfect fluid with energy density 
p and pressure p. Though a constant density fluid is not realistic, it represents a benchmark 
and in GR a great deal of interesting information on the possible stable gravitational bounded 
configuration can be derived. Different from GR, because of the presence of the Goldstone fields 
even for a constant density fluid the resulting Einstein equations are difficult to solve analytically 
(the full system of equations in presence of matter is given in Appendix [Cj and we have to rely 
on perturbation theory. Moreover, even the linearization procedure is not straightforward and one 
has to set up the nondemocratic weak-field expansion as discussed in the last section. 

The starting point is the expansion in the weak-field e parameter 



J = l + eJ (1) + --- , cf> = br + ecf) {l) + ■■■ , p = e p + • • • 

K = l + eif (1) + ••• , h = e 1/2 + • • • , p = e 2 p {1) + 



(13) 



where po is the constant density of the fluid. As discussed above, the different parametric size 
of h is necessary to have a consistent expansion of the t-r component of Einstein equations. Diff 
invariance of the matter action alone leads to the conservation of the matter EMT. Then, by 
expanding the metric, the pressure is of order e 2 and can be neglected at leading order. 

The linearized solution contains four integration constants, two are set to zero by imposing 
regularity in r = 0; the last two, together with S and M are determined by matching the interior 
and the exterior solution at the star radius R. For brevity, we only give here the solution for the 
case A = and 7 < 0. The final expressions are 



7(1) _ GM Q 

J - R 



" 1 + <S h ' I 2 -1-2 " 5( 7 -3)( 7 -4) + ( 7 + 3)(2 7 -l)(7 + l)(7-2) 



R l? ; ' V °V(7 + l)(7 + 3) 2(7 + 3) 2(7 + 1)7' 1 > 

where Mq = | ir R 3 po is the bare mass of the star, £ = r/R and finally A = JK so that A^ 1 ) = 
jM+id 1 ). 5 

The parameter p is a new important mass combination defined in terms of the masses around 



flat space, (12): 



3 m\ — to 2 , (m^ — 3 m 2 " 1 
m\ — tuq (m| — m 2 .) 



5 For 7 = — 1, —3, eq. ( 141 becomes singular and these cases need a special treatment given in appendix |f| 



s 



Note that in general fi 2 can be positive or negative. The function h always appears as h' 2 and is 
obtained in terms of d>^~' and A^: 



h 
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12 m 2 r 2 



(A 



(i) 



r j(D'_ j(D)_ j(D + (^) (i) 



(16) 



The last unknown, the pressure p(r), is found by using the matter EMT conservation, and is 
expressed in terms of J and an integration constant po, fixed uniquely when defining the radius R 
as the point where p(R) = 0. The pressure is of order e 2 but the ratio ^ is of order one, therefore 
we have 



p{r) 
Po 



J(R) 
J(r) 



nl/2 



2 L 



J (1) (i?) - J (1) (r)l +0(e) 



(17) 



By matching J, J' at r = R with the exterior vacuum solution (11) (with A = 0), we find the two 



exterior integration constants M and S in terms of the parameters of the star: 

8u 2 R 2 1 _ 2Au 2 M R 1 ^ 



M = M, 



1 



5( 7 + l)(7-2) 



(18) 



( 7 -4)( 7 + l)(2 7 -l)( 7 

We thus find that the star acts as a source for the new term, 5^0, and that the bare mass Mq 
is renormalized. For 7 < — 1, both S and the mass shift AM = M — Mq have the same sign of 
fi 2 . Thus, for u 2 > both corrections are positive, M > Mq and S > 0, while for fi 2 < both the 
corrections are negative and then M < Mq and S < 0. 

The deviation from GR is measured by S which scales as i? 4-7 ; as a result, bigger self-gravitating 
objects produce larger deviations. The difference between the gravitational mass seen by distant 
observers and the bare mass, AM = M — Mq, can be traced back to the Goldstones' energy density: 
we have, using the equation of motion (EOM), 

S(7 + l) 



,7-2 



4tt 



J(r) 



(19) 



Thus in the exterior region the energy density of the Goldstones is positive (for fj? > and in the 
range 7 < —1 we are considering). Moreover, having a regular solution in all the spacetime, the 
Komar energy can be computed as an integral of the total energy momentum tensor (matter + 
Goldstones) over a t = const 3-ball of radius f: 



d 3 x 



rpV 



2 M 



/ r \ 



en v = Mq + AM + Ci m z R° Po {-) 



7+1 



' t=const 

with C\ an irrelevant constant. For 7 < — 1, the energy is finite and equal to M in the limit f 



(20) 



00. 



Since via fi 2 both S and AM are proportional to the graviton mass scale m 2 , we conclude that 
for m — > the deviations from GR disappear, both in the interior linearized solutions and in the 
exterior exact one: 



J« -+ J, 



(i) 
GR ' 



J -> ^GR • 



M -> M , 5^0, 



(r < R) 
(r > R) 



(21) 
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fi~<0 dashed green 
fj~>Q dotted red 
GR solid black 





GR solid black 




/i">0 doited red 




f.r <0 dashed green 
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FIG. 1: Plot, in the linear regime, of J — 1 and p/po in the case of massive gravity (14 1 and (17), both for 
\j? positive and negative; also the standard GR case is shown for comparison. The numerical values used 
are : y?R 2 = ±3.06, R s /R = 8.6 x 1CT 7 and 7 = -2.5. The units are such that GM /R = 1. 



where GR indicates the corresponding expressions in general relativity. In the limit of vanishing 
graviton mass, GR is smoothly recovered even in the presence of the Goldstone fields (f>, h that 
do not got to zero. Thus, there is no discontinuity for m — > in the spherical star solution, 
and there is no sign of the Vainshtein scale ry that in the Pauli-Fierz theory invalidates standard 
perturbation theory at scales r < ry. Here, the standard weak- field expansion is never valid, while 
the nondemocratic expansion is valid at all distances; in other words, ry = 00. 

The above conclusions apply to realistic weak-field stars like the Sun. For the strong-field regime 
the inner star solution can be found numerically, and we discuss below some of its properties. Of 
course the present linearized solution coincides with the exact numerical one when the density is 
sufficiently low. In Fig. [T] the GR expression of J and K are compared with the ones found here 
for massive gravity. 



C. Internal pressure and a bound on the graviton mass 

The dimensionless parameter p?R 2 plays a key role in determining the size of the deviations 
from GR. In particular it is interesting to investigate both the behavior of the internal pressure 
and of the total mass M, because their renormalization is sensitive to p?R 2 and may turn negative 
when fi 2 R 2 is order one. 

Clearly the stability of the matter system requires the pressure to be positive. One can see that 
it is enough to impose this condition at two regions, near the center and near the surface of the 
star. In particular, at the center of the star the requirement of positive pressure is 



p(0) 
Po 



GM 
2R 



1 



16/i 2 R 2 (11 - 2 7 ) 
5(2 7 -l)( 7 -4)( 7 -2) 



> 0. 



(22) 



while to have positive pressure near the surface, where p(R) = 0, its derivative has to be negative: 



p'(R) GM 



Po 



2R 



16n 2 R 2 (7 + 2) 



5(2 7 -l)(7-2)( 7 -4) 



< 0. 



(23) 
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FIG. 2: Excluded regions, in the (n 2 R 2 -j) plane, from positivity of the pressure [p(0) solid line, p'(R) 
dashed line] and of the renormalized mass M (dotted line). The central white region is allowed. 

On the other hand, the requirement that the renormalized gravitational mass should be positive is 

o 2r>2 

M = M 1- , . /" >0. (24) 

°[ 5( 7 + l)(7-2)J 1 ' 

Thus, when 7 < — 1, we have Sp(0) oc fj, , 6p'(R) oc [i 2 sign(7 + 2) and SM oc — fx . As a result, for 
positive /j, 2 the bound on fj? comes from p'(R) (for 7 < —2) and M, while for negative /j, 2 the bound 
comes from p(0). This covers the whole range of fi 2 . In Fig. [2] we show the range of parameters 
spanned by such a requirements. The stronger bounds are due essentially to M and p(0). 

Assuming no cancellations, or in other words that all mass scales in the potential are of the 
same order, the limit on /i 2 can be interpreted as a limit on the graviton mass scale m 2 . 

The typical limit is therefore m 2 ~ \fi 2 \ < 0(1)/R 2 ; therefore, real stars with e.g. R ~ 10 5 km, 
require m < xlO~ n eV. Considering the Sun, for which R ~ 7 x 10 5 Km, and for which the 
central pressure may not deviate more than few percents from the standard value [21] , we have 
m < xlO _13 eV. However, since the limit does not depend on the mass of the body but only on 
its radius, it can be applied also in the extreme weak- field limit (R 3> R s , see next section), for 
instance to very large and rarefied objects. In fact, larger objects give rise to stronger limits; for 
instance, the stability of the largest bound states in the Universe, R ~ 1 -j- lOMpc, gives a limit 
m < 10- 28 - 29 eV. 

Though the above considerations are strictly valid for a spherical symmetric and constant density 
body, we believe that for more realistic configuration the above limit can capture the correct order 
of magnitude. For other independent limits on graviton mass based on pulsar timing and CMB 
polarization effects see [19j. 6 



6 One should note that the usually quoted limits on the graviton mass, like some reported in the PDG, should 
be taken with some care. For instance, the strongest reported limit [50], m < 10 -32 eV, is derived assuming a 
Yukawa-like potential. Therefore, in the class of theories we are considering, such a limit does not apply. 
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Let us also comment that the positivity of the renormalized mass M is not strictly required. 
In fact, a negative M (and 7 < —1) would lead to a repulsive gravity at large r (only) while 
the body would still be in equilibrium, having positive pressure. Such and exotic behavior of the 
matter-Goldstones system may deserve further attention. 

Regarding the validity of the bounds described above, an important discussion is in order. 
Clearly, the bounds hold, assuming that the description in terms of gravity coupled to a fluid 
is valid at the scales of interest, e.g. galaxies or a cluster of galaxies. This issue is far from 
being trivial, since, e.g. in standard gravity, the application of the same theory at different scales 
is guaranteed only in the Newtonian limit by the linearity of the field equations. In nonlinear 
cases, like for strong field or in the cosmological evolution, the description in terms of averaged 
quantities is expected to need effective corrections, which are under active discussion (see e.g. |22j). 
In the present theory, two comments can be made. First, the description of matter in terms of a 
gravitationally bound fluid is appropriate, because the typical interparticle distance is small with 
respect to the size of the body and to the range of the gravitational force, which we recall is always 



infmite[(the Newtonian term in Eq. (11)]. On the other hand, the averaging of the gravitational 



field equations is non trivial, because, as it is shown by Eq. (16), one of the Goldstone fields 
enters quadratically, (h 1 ) 2 . (accordingly, h! is of order yfe ~ 1/\A" h 1 the exact and semilinearized 
approach that we described). Because of this nonlinearity, one expects a violation of superposition 
of multiple solutions, even in weak- field regime 7 Notice that this nonsuperposition has to arise 
also in the Lorentz-invariant version of massive gravity [6], where one of the fields, essentially h' , 
enters quadratically and has a 1/ y/r falloff . We expect in fact this to be a generic phenomenon in 
massive gravity. What we can conclude is that if one assumes the present description of gravity 
to be applicable at some scale of interest, i.e. galaxy clusters, then stability of the system leads to 
the strong bound derived above. 



IV. PERTURBATION THEORY AND BEYOND 
A. Validity of Perturbation Theory 

Let us now study the validity of perturbation theory in both the inner and outer region of the -star'- 
. In GR it is well known that perturbation theory can be used when r/R s 3> 1, with R s = 2GMq 



the Schwarzschild radius. By inspection of the expression for the metric perturbation ( 14 ) inside 
the body it is clear that perturbation theory is valid when 

R S <.R, fi 2 R s R<.l. (25) 

In this section we use fi 2 as \[i 2 \ always positive defined. 



7 Not to be confused with the violation of Birkoff theorem or of Gauss law, which can be violated even in linear 
equations, when departing from purely Newtonian behavior. 
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If we were unable to find an exact solution in the outer region we could have set up the very 



same nondemocratic perturbation scheme as in (13). In fact, solving the linearized equation we 
find 



Ji 



+ /ir 2 +j2r 7 + /2r i-7 ) K (i) = _jW j 0« = / 1 r + / 2 r- 



7 • 



(26) 



where j\ 2 and /1 2 are integration constants. When 7 < —1, to be in the weak-field regime at 
large r we have to set f\ = fz — 0, and we get the linearized version of the exterior solution. The 



expansion is valid when jx/r <C 1 and *C 1. Using the values (18) of jx and j 2 obtained from 

(27) 



the matching with the interior solution, we get, for r > R, the conditions 



The first condition makes sure that we are away from the would-be horizon at tq and it is trivially 
satisfied in the presence of the interior part of the solution with R > R s . Notice that in the 
range 7 < — 1, thanks to (25), also the second condition is automatically satisfied. Then, once 



perturbation theory is valid in the interior of the body it can also be used in the exterior part. 



When (25) is not satisfied, the determination of S and M cannot be done using perturbation 



theory. In this case one can solve numerically the equations in the interior part and match the 
numerical solution with analytical exact solution that we have found in vacuum. 



It is worth stressing that even within the perturbative region (25) sizable corrections to the 



gravitational mass M are possible. This is the case for large bodies with R of order 0(1)/ fj, . 



In this case, the second relation of (25) becomes equivalent to the first one R s <C R, but the 
dimensionless quantity fi 2 R 2 appearing in M can be sizable. 



B. Beyond Perturbation Theory 



When the equations in (25) are violated, perturbation theory cannot be trusted anymore and 
one needs a different tool to see what happens in this regime. We have solved numerically the 
Einstein equations in the interior, still modeling matter as perfect incompressible fluid. Figure [3] 
shows M and S computed numerically for different values of the graviton mass and matter density. 
From the numerical analysis it is clear that the linearized matching captures the basic features of 
the geometry. Indeed, the difference with the linear predictions are rather small and below 15 % 
for R s /R ~ 0.5. 

The numerical analysis confirms the absence of a discontinuity for vanishing graviton mass, as 



discussed above in the perturbative analysis [see ( 21 )] . Lowering the values of m, while keeping fixed 



the density, we have found a behavior compatible with S, AM — > 0, up to the moderate strong-field 
regime R s /R ~ 0.3, and down to small graviton mass scale mR ~ 10 -4 . The numerical analysis 
thus shows no sign of discontinuity in m. 

The remarkable feature that larger bodies gravitate more survives at the nonperturbative level. 
The relative mass renormalization for a body of size R ~ mT 1 is large, AM/M ~ 1. In Fig. [4] 
we extend the perturbative bound on the dimensionless combination mR to the region where 
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FIG. 3: Numerical results for AM = M — Mo (left) and S (right) for different values of the matter density 
and of the graviton mass scale /j, 2 R 2 = 0.1,..., 1 (lower to upper). The dashed lines correspond to the 



perturbative values; ( 18 1 



R s /R ~ 0.5. The result is that toward a strong field, the Goldstone pressure becomes even more 
negative, pushing the limit on mR even lower. Of course real heavy objects are always small, 



therefore the bound given above in Sec. IIIC for large weak-field objects is more stringent. 



0.5 , 



0.4 - 



0.3 



0.2 - 



0.1 - 



0.05 0.10 
2GM/R 



FIG. 4: Numerical exclusion plot of the graviton mass [from p(0) > 0] toward the strong field. Here 7 is 
fixed and the perturbative limit from (22 1 is fi 2 R 2 < 0.338 (dashed line). 



Finally, the fact that the central pressure receives a negative contribution from the Goldstone 
EMT; points to the possibility of relaxing the Chandrasekar bounds on the density which, we 
recall, arise in the strong-field regime when the central pressure diverges. There may be then a 
regime where the negative pressure of the Goldstones balances the high gravitational force, allowing 
for compact heavy objects with arbitrarily high mass. We could not reach this regime with our 
numerical analysis that becomes arduously difficult near the origin, when going beyond R s /R ~ 0.5. 
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V. MASSIVE GRAVITY BLACK HOLES WITH HAIR 

It is known that in GR the gravitational hydrostatic equilibrium cannot be attained for arbitrary 
M and R, for instance a body of uniform density cannot have a mass larger than M = 4R/9, and 
for a fixed radius, this limit applies to any equation of state. The limit is crucial, because when 
a star in its late evolutionary state has exceeded it, the gravitational collapse cannot be stopped 
and a black hole is expected to form. 

In GR, the spherically symmetric (uncharged) black holes are parametrized by its mass M. In 
massive gravity, as we showed, they are parametrized by an additional parameter, S. Suppose then 
that the late stage of the evolution of a star is described by our exterior metric 

ds 2 = -Jdt 2 + Kdr 2 + r 2 dn 2 , J = 1- 2G E + 2GSr^ = K~ 1 . (28) 

r 

and for definiteness let us consider here the case with M > 0. In general, the hyper-surface at 
r = rjj where J(th) = is a Killing horizon: at r = rjj the norm of the timelike Killing vector 
d/dt vanishes. If J is never vanishing the singularity at r = is naked. This happens when S is 
large and positive, S > -^q-- When S is negative J has a single zero. Finally, there is a range of S 
where J has two zeroes or one, namely < S < -^r. When two zeros of J are present, typically 
the outer Killing horizon is an event horizon for the black hole. For instance, in the case 7 = — 2 
the solution is formally the same of a Reisser-Nordstrom black hole. To disentangle the causal 
structure in the general case one has to try to maximally extend the solution found using a set of 
Schwarzschild-like coordinates that fail at the outer Killing horizon. We leave the detailed study 
for a future work. 



An other evident feature of the black- hole solution (28) is that it is one of the few examples 
of exact black holes with hairs in four spacetime dimensions. Their existence is considered to be 
forbidden in various classes of theories, by "no-hair" theorems (starting from |23j). These are 
usually based on the positivity of suitable volume integrals, and often rely on linearity of the 
field equations. For instance, a no-hair theorem for the Fierz-Pauli massive gravity was described 
already in |24j . where only the quadratic FP mass terms were included. As we have proven here, 
the full theory of nonlinear ly interacting massive gravity on the contrary generates hairs, extending 
in the asymptotic region. These hairs, depending on the value of 7, can even dominate over the 
Newtonian term, and even make the solution not asymptotically flat (7 > 2), leading to a type of 
furry black hole. 

The situation is thus similar to known cases of black holes with hairs in non-Abelian gauge 
theories (colored black holes). However, classical non-Abelian fields are actually not observable 
(due to confinement). The example provided here would be the first example of classical hairs 
from nonlinear ly interacting fields. 

In the literature there are a number of no-hair theorems which also consider specifically the 
presence of extra scalar fields (see for instance |28j ) that, in principle apply to our case. However, 
"technical" simplifying assumptions on the form of scalar field Lagrangian are often made, rendering 
the above results not directly relevant to our nonminimally coupled and interacting set of scalars. 
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We stress that the fields associated with these hairs (both the metric and the Goldstones) are 
nonsingular in all the space, except at the origin, where in any case the singularity is expected. 
Also, in the smooth star solution that we constructed, these fields are never singular. 

Clearly, whether hairs will actually be present in real black holes, i.e. whether S will be nonzero 
at the final stage of collapse, will depend on the collapse history, which is surely worth of a separate 
study. For instance, the fact that in the star solution the corrections to GR are proportional to 
the star's surface area, i.e. AM, S ~ R 2 (but only in weak-field) may lead one to think that 
S — > during collapse. Or on the contrary one may think that the collapse itself may be stopped 
or reversed leading to a final "remnant" with a non-zero S. We stress in fact that the Birkhoff 
theorem is violated in this theory, and therefore the real dynamics is only constrained by the total 
energy (not just the Komar integral). A related important issue is whether the exact solution is 
stable under perturbations: if on one hand the non positivity of the Goldstone potential T hints 
at instability, on the other hand we recall that at least in perturbation theory, due to Lorentz 
breaking, there is no badly behaving propagating mode. We leave the study of stability of the 
black hole and of the possible scenarios of its time evolution for a separate analysis. 

Finally, the classical hairs, as realized in our solution, seem to probe the structure of the black 
hole inside the horizon, and thus may alter already at classical level its thermodynamical properties. 
The power-law hair SV 7 , extending toward the singularity, may also be sensible to quantum gravity 
effects, and thus be a probe of the scales of UV completion. 8 

VI. CONCLUSIONS 

One of the difficulties of massive gravity theories is that perturbation theory can be very tricky, 
making hard to extract solid phenomenological predictions. This is why it is crucial to find exact 
solutions at least in highly symmetric configurations. The extension of the familiar Schwarzschild 
solution is a step forward for testing massive gravity. 

In this work we have addressed this problem in a wide class of promising massive gravity theories, 
where the only propagating degree of freedom is a single massive graviton [H] . In the Stiickelberg 
spirit, general covariance is restored by a set of suitable (Goldstone) fields, which are nonlinearly 
interacting with the metric field. 

We have determined an exact class of black-hole solutions describing the exterior gravitational 
field of spherically symmetric compact bodies, and described its matching with an interior part, 
describing the structure of a self-gravitating body. This solution was found analytically by a 
suitable nondemocratic weak-field expansion. 

For the exterior part of the solution, we have found that gu = g^r as m GR. Incidentally this 
implies that the PPN parameter, which measures the difference between the gravitational potential 



Such considerations have also been put forward in |27| . in connection with the Lorentz breaking. In that case 
however, the focus is on states propagating with different speed of light, thus probing the interior of the horizon; 
here, we remind that there are no other propagating states beyond the ordinary graviton. See also |25j in the case 
of the Lorentz breaking ghost condensate, and for consideration of LB theories with instantaneous fields. 
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felt by a massless and a massive particles, vanishes identically, and as a result the light-bending in 
LB massive gravity coincides with the one of GR. 

The differences with respect to GR appears inside gu, where in addition to the Schwarzschild 
term, a new nonanalytic one appears, gu = 1 — 2GM/r + 2G5r 7 , where 7 is a constant that 
depends on the coupling constants of the theory and S is a new free integration constant. The new 
term Sr 1 represents a genuine black-hole hair due to the interacting Goldstone fields. Depending 
on the value of S there are a number of possible scenarios ranging from a naked singularity at the 
origin (S large and positive), to a strong gravity at large distance (for 7 > 0), or to a standard 
asymptotically flat behavior with a finite Komar mass equal to M, when 7 < — 1. We stress that 
the phase of Lorentz broken massive gravity we are considering, the range of the static gravitational 
potential is infinite even in the presence of massive deformation. 

For a realistic spherical body like a star, the Schwarzschild mass M and the new parameter S 
can be computed (in the physical case 7 < — 1) by matching the exterior (exact) solution with the 
interior (weak-field) one. This leads to two results: 

• The usual "bare" mass of the body Mq is renormalized, M = Mq(1 — cq ^ 2 R 2 ), where R is 
the radius of the star and cq is some numerical factor. 

• S is turned on by the presence of matter, S ~ Mq /i 2 R 1 ^ 1 ' . 

Here /j 2 is a combination of graviton mass scales, which in principle may be negative. The sign of 
fj? controls to what extent the -bare- mass gravitates as seen by a distant observer: when /j, 2 < 0, 
M is larger than the 'bare' mass; on the other hand when fi 2 > the body "degravitates" . The fact 
that the gravitational field depends on the shape (size) of the body is a signal of the violation of the 
Gauss law, produced by the nonlinear interacting Goldstone fields. In both cases it is striking that 
AM = M — Mq is proportional to the body surface R 2 ; thus, larger bodies (de-)gravitate more. 
Surprisingly enough, degravitation can be so large that M < while the body internal structure 
is still in equilibrium. 

Also the internal pressure of the body is subject to a similar renormalization. By applying this 
theory of gravitation to known self-gravitating structures, one can then derive a bound on [i 2 R 2 
by requiring positivity of the pressure, or in other words by requiring their equilibrium. Normal 
stars like the Sun pose very loose constraints on /j, 2 , but the effect is stronger for larger objects, 
so from the largest (and less dense) bound states in the Universe one may pose a strong limit of 
order fi < io~ 28H ~ 29 eV. If the various Lorentz-breaking masses are to be of the same order, this 
translates into a strong constraint on the overall graviton mass scale. 

The renormalization of mass and pressure, as well as the hair S r 7 are directly produced by the 
presence of the fluid made of Goldstone fields. However, the Goldstone fluid is seeded only by the 
presence of matter (at least in weak-field), and the effects disappear for Mq — > 0. In other words, 
there is no smooth spherical body made of Goldstone fluid only. Similarly, the deviations from 
GR disappear for vanishing graviton mass scales, m 2 ~ fi 2 — > 0, showing no sign of discontinuity. 
This result is exact in the (nondemocratic) weak field limit, but a numerical investigation of the 
solutions in strong- field conditions was performed. Up to 2GM/R ~ 0.5, the analysis confirmed 
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the absence of discontinuity both in the limit of vanishing matter or vanishing graviton masses. 

We did not address the difficult problem of the stability of the proposed solutions, which is 
very interesting in view of possible phenomenological applications, and also in view of the current 
conjecture of the instability of configurations with hairs |29| . Notice that because the Birkhoff 
theorem does not hold for the gravity modification we have considered, the study of stability is 
rather complicated. 

Finally, let us also emphasize that an analysis of the propagating modes beyond the linearized 
order is still missing and should be addressed for these theories to be considered viable. We leave 
these important studies for further work. 

Acknowledgments. During this work D.C. was partially supported by the EU Contract No. FP6 
and the Marie Curie Research and Training Network -UniverseNet- (MRTN-CT-2006-035863). 

Note Added. The authors of [TB] consider the same problem of finding a spherically symmetric 
solution in Lorentz-breaking massive gravity for a particular potential. According to that paper, S 
is always zero and <j) = r for r < R. However, this case corresponds to a body made of cosmological 



constant [see Eq. (C3) in Appendix [CJ for which A = 1 even in the interior, and <j> = r everywhere. 
Of course, this is not true for a generic and realistic kind of matter, and S can be different from 
zero as shown in this paper. 



Appendix A: Various approaches to Massive Gravity 

Building masslike terms for the graviton requires the possibility of scalar combinations of metric. 
This can done in an elegant way by introducing a metric Gab in a fictitious manifold At, the set 
of four "scalars" can be interpreted as a mapping of the physical M space into the fictitious one, 
<3? : Af — > A4, in terms of coordinates — > $ A (x). The metric in the fictitious space then can 
pushed back to M rendering available the basic tool for constructing diff invariants generalized 
mass terms 

d$ A d$ B „ , . 1 , 

92 ^ = ^^ Gab - (A1) 

The new metric gi transforms as a standard tensor and can be used to construct nonderivative 
interaction terms by introducing 



XZ=g 2 ^g av , r n = Tv(X n ) = (XX- (A2) 
A typical interaction term will be of the form 9 

S int = m 2 M 2 P j d A x^gV(n, T 2 , 73, n) (A3) 



One can also consider combinations involving y n = Tr(Y n ), Y = X 1 . 
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The metric Gab may or may not be considered as dynamical variable. When G is dynamical, the 
full action is invariant under two separate diffs Diff 1 x Diff 2 corresponding to the physical space 
Ai and to the space Ai. If G is frozen to some background value, the invariance is broken down to 
the set of diagonal cliff Diff d . The pointwise identification of the two manifold is obtained imposing 
$> A (x) = 6^, the so called unitary gauge (UG). In the UG, diff invariance is broken down to Diff d 
when G is dynamical or completely broken when it is a frozen background. As an example, taking 
V = — T2 — 6ti + Co, in the UG, setting g^ u = + h^ v and Gab = Vab, expanding up to the 
quadratic order one gets the Pauli-Fierz model. 

Appendix B: Goldstone Energy Momentum tensor 



The Goldstone energy momentum tensor for a generic F{X , W l i) is given by 



T g = m 2 M 2 P l Fg» v + 2F X + 2F t 



X -2 yi yj yOO _ #-1 (yiyJO^ + yjyiO^ + yi^ 



)J,V 



with Y AB = d^ A d v $ B and 



dF 



F x — o ) 



dF 



(Bl) 



(B2) 



Since Fij a 3 x 3 matrix that depends on W 13 , it can always be written as 

Fij = F 1 & + F 2 + F 3 {W 2 f , (B3) 
with F\ , Fi , F3 as the scalar coefficients that depend on the explicit form of F. 



Appendix C: Einstein equations 

We give here the full set of Einstein and conservation equations. 

The matter EMT conservation gives a first order differential equation for the pressure 10 

p'+ J ^ = o. (CI) 

The off-diagonal component tr, simplifies considerably for the class of potentials considered and, 
when h! 7^ and gives an equation that can be solved for (f)' 2 : 

, 2 = r 4 6 6 A(2r 4 6 4 /3 2 -2r 2 6 2 /3 1( A 2 + /jo0 4 ) 
^ r 4 6 4 ai0 4 -4r 2 6 2 a 2 6 + 6a 3 </> 8 1 ' 



10 The same equation is also valid in GR, where in addition also J'/ J can be eliminated using the rr Einstein equation 
leading to the Oppenheimer-Volkoff equation, which generalizes the equation of hydrostatic equilibrium. 
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Even independently from the potential, using the tr equations it turns out that the difference 
between the tt and rr components of Einstein equations becomes 



ZV _ r (p + p) 
A 2 ~ 2M| J ' 



(C3) 



We note that (C3) holds also in GR. It is interesting to point out that in the presence of matter 
in general A = JK 7^ 1 , unless p + p = that corresponds to the vacuum or to a body made of 
cosmological constant. When p + p = 0, like in the exterior part of the solution, = br and the 
Goldstones' EMT is rather simple: 



2 (ai - 4a 2 + 6a 3 ) m 2 M 2 a 2 J - a 2 + J 2 (tif 



T/ = f{r)5t 



f{r) 



a 2 - J 2 (/>') 



(C4) 



The tt component of the Einstein equations depends on the particular structure of the potential. 
Its general form can be given as 



1 (J' J\ p 1 m 2 A 
+ ~ 1 + 2.U 2 + T 2 + 



A V r 



0. 



(C5) 



where H and L depend on T . Given the spherically symmetric ansatz ([5]), the explicit expression 
for X is (notice that X cx m 2 ) 

K-Jh' 2 



X 



A 



and one can solve rtC5| for {tif. Explicitly for our potentials 



(h 



" 2 -= a 2 J- 2 



1 



Jm r (a± — 4a 2 + 6(^3) 



r (J'(r) + m 2 r (a\ — 4a 2 + 603)) + J — 1 
Finally, the 96 component of Einstein equations can be written as 



J" + J' 



, (2 _ r A (p + p) 
AJMp 



12m 2 A 



di{4>) 
X 



+ d 2 



A(3p + p) 
2M 2 



(C6) 



(C7) 



(C8) 



where again the functions d\ and d 2 depend the choice of J-. 



Appendix D: Behavior of the solutions in the limit m — » 



An anomalous behavior of the solutions in the limit m — > can take place due to the form of some 
of the EOM's 



m 2 f{x,y) = 0, 



or 



x = m 2 g(y) 

y = irK x ,y) 



(Dl) 



where /, g and k are smooth functions of their arguments, so that the solutions obtained for 
m = are in general different from the ones obtained by taking the limit m — > of the generic 
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solutions. Equation T gtr = is precisely of the first form, while the tt and rr components of 
Einstein equations are of the second form. For the Ssake of compactness in this section all the 
explicit formulas refer to the special case: ao = 2, (3o = — 1, /3i>o = 0, 02 = «4 = 0. Assuming 
m / 0, we can solve the equation T gtr = for 

2r 3 ( 7 - 2) 



J" + J' 



r 4 ( 7 - 2) - ~/<p 4 



r A (p + pp) 
4JM 2 



+ 



and Eq. (C5) for X; then (C8) becomes 
2J 7 c/> 4 



r 6 (7 — 2) — r 2 j(p 4 



12m 2 (7-2)A 



r 6 ( 7 -2) -r 2 7</> 4 



A [ 7 4 (r 2 (p + gg) - 4M 2 ) - r 6 ( 7 - 2) (3p + go)] 
2M 2 r 2 (r 4 ( 7 -2) - 7 </> 4 ) 



(D2) 



0. 



The GR equations (obtained imposing m = from the beginning) are instead given by eqs (CI) 

'2 rA{p + p )\ A(3p + p ) 



and (C3), and 



J" + J' ( - 

r 



AJMp 



2M 2 







In vacuum, the difference between Eqs. (D2) and (D3) is not zero even if we put m 

2 7 4 (rJ' + J-A) 



in (D2): 



(D2)-(D3)\ i 



in vacu um 



r 4 (7 



7 </> 4 



(D3) 
directly 

(D4) 



As a result, the equations themselves are discontinuous and we expect a discontinuity in the space 
of solutions. 

Incidentally, as it is easy to verify, the standard Schwarzschild solution of GR, i.e. J = 1 — 
2GMq/t and A = 1, satisfies exactly the above expression. This means that the Schwarzschild 
solution is also a solution of the Massive Gravity equations, but these can have new solutions, 
corresponding to a nonvanishing Goldstone EMT. This is indeed clear in the exact solution given 
in the text, where a new term with a new integration constant; SV 7 ; is present. The discontinuity 
in m can also be understood from the exponent 7 which is given by a ratio of mass parameters of 
the Goldstone potentials, and as such persists in the m — > limit. 

The question is then whether for realistic star solutions the constant 5" is vanishing or not for 
m — > 0. In the text we have shown that in the weak-field regime the solution is smooth in m: both 
AM, S — > so that the solution reduces to the standard GR one. 

In general, inside a medium, the structure of the EOM is the following 



S-kG 



(D5) 



The presence of the Goldstones' EMT introduces more equations than in GR. For instance, in GR 
both the Einstein tensor and the matter energy momentum tensor are both diagonal; it is not so 
when the Goldstones fields are introduced: 

6Jti (2r A JK+ (r 4 ( 7 - 2) - 7< /> 4 ) 



k'2 



rp t 

1 9r 



-m 



) 



(D6) 



r 4 j(K-Jh' 2 ) 2 

The Einstein equation tr leads to T g t = that allows us to solve for <j>' (Of course we assume that 
neither J nor h! are vanishing) 

,,2 2r 4 JK 



r 4 ( 7 - 2) 



(D7) 
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In GR, introducing the function v s { r ) by 

Kir) = 



2G_ 



-i 



(D8) 



the tt equation can be solved in terms of the integral of the matter energy density. In our case the 
tt equation is more complicated, 



2 3r m 
vnr p H — 



(7-2)<^ 2 
r 2 7 



2rJ 



7J(r - 2Gu s )k 



/2 



1 



r7 



(D9) 



Notice that again the previous equation has a continuous limit when m — > and reduces to the 
one of GR unless J, h' or K are singular when m — > 0. The discontinuity turns up when one uses 
the 86 equation to eliminate h': 

m 2 U n 



h 



)2 



K 

J m 2 T-Ldi + Ud2 



(D10) 



where Hi are suitable functions of J, K and (ft that do not depend explicitly on m (they will 
depend in general on m implicitly). Replacing h' 2 in (D9) by using |d|, one finds 



a 2 2,[2r 2 (U dl m 2 +U d2 )K + U n m 2 {r 2 1 -{ 1 -2) ( t ) 2 )] 

GHnl 



(Dll) 



As a result, when one takes the limit m — > 0, supposing the J, K and are regular, Eq. (Dll ) differs 
from GR if lim m ^o^d2 7^ 0. At the linearized level, using (14) one finds that lim m _ > o^d2 = 0, 
explaining why the solution is smooth in m. 

The question whether the strong-field solutions are continuous for m — > can at this stage only 
be addressed numerically. Up to the field intensity of R s /R ~ 0.5 we could verify that the solution 
is actually continuous. 



Appendix E: Gravitational Energy 

The gravitation energy can be evaluated by the standard Komar mass |17j . In the presence of a 
timelike Killing vector £^ the gravitational energy is given by 

E = --±-( ^-pd 2 xn^V^u , (El) 

where (?Ej, with unit normal v^, is the boundary of the spacelike 3-surface t = const with unit 
normal n M , finally p^ v is the induced metric in <9£t. In our case we take <9£t to be the 2-sphere 
t = const: r = f of large radius r and the Killing vector is £ = d/dt. We find 

E f = M + S~if 1+1 . (E2) 

When 7 < — 1 we can take the limit f — > oo and -Boo = M. When 7 > — 1 the Komar energy is 
infinite. A detailed study of the Hamiltonian approach and in particular of the case with 7 > — 1 
will be given elsewhere |18j . 
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Appendix F: Special cases 



1. The case 7 = — 1 



In the case with 7 = — 1 we have to change also the exact vacuum solution for the J function: 



_ , 2GM 2GS, (t 

J = l + log „ 

r r \R 



(Fl) 



the solution of the linearized equations inside the star of constant density po read 
j(i) GM 



R 



A(1) = 3GMq ^ 2 _ ^ 

* (1 > = ^^fe S -f (1 + 2 log*)] 



(F2) 



where £ = r/R and jt*~ = _i = (— ai + 16a2 — 42a3 + 12/32 — 36/?3) "i 2 - The matching at the boundary 
r = R gives the parameters M and 5 in terms of the star bare mass and radius 

13 f?2-2 \ q 8 

lOO^^-V' S = ~5 
The (Komar) energy inside a large shell of radius f is given by 



M = M 1 + 



(F3) 



& = Mi 



1 



12 



2b R 2 ^ = _ l +~-R^i 1= _ l \og yR 
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j2,.2 



(F4) 



and diverges as a log in the limit r —> 00. 



2. The case 7 = 3 



The linearized solution given in the text (13) is modified for 7 = —3. The exterior solution has 



of course the same form, while the interior solution is drastically different: 



GM 

R 
3GM ( 



r - 3 + i? z A 



7=-3 



f - or + -^riog£ 



V 1960^ 8 
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/>2 



(e 2 - 1) 



(F5) 



(i) = 36G^ [e + f8( _ 1 + 21ogO]) 



where £ = r/R and \xt f= _^ = (— 9qi + 56a2 — 114«3 + 20/^2 — 6O/33) m 2 . The energy is finite and 
the parameters of the external solution are given by 



.\f=M [l + ^R 2 fiZ = _ 3j 
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S = -^—M R 4 j2 2 



7 =-3 



(F6) 
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